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SUMMARY 


The present methods available for the determination of the 
magnitudes of the vibrational relaxation times of molecules are 
discussed and the recently developed theory of Schwartz, Slawaky, 
and Herzfeld is used to compute the variation with temperature of 
the collisional vibrational excitation probabilities of mixtures of 
N> and O>» and of Eo and HE. The results are used to estimate 
the extent of vibrational temperature lag in the hydrogen-fiuorine 


rocket motor. 
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l. INTRODUCTION 


The internal energy, E, of a gas may be written, approxi- 
mately, as the sum of energies of translation, rotation, vibration, 


and electronic motion. That is to say, we write that 
E = Ey anelT) + Eve CT.) + En (T) +E (TD, a) 


where = a Ty» and Te. represent, respectively, the absolute temp- 
eratures associated with the translational, rotational, vibrational, 
' and electronic energy states. Under some conditions these energy 
states will be in equilibrium such that T=T)>T,- e° If, in the flow 
of a gas, ternperature changes occur at rates sufficiently fast to ex- 
ceed the time required for these internal energy states to readjust 
themselves, the gas will depart from its equilibrium partition of 
energy. 

The term vibrational relaxation time" has. come into use as 
@® measure of the time required to establish equilibrium with respect 
to the internal vibrations of the molecules. The other internal states 
adjust themselves so rapidly that we can usually neglect their relaxa- 
tion times. The vibrational relaxation time, T , is defined by the 


relation 


d Eyip(T//at = v[E a (T) 3 E 4 (T.)| : (2) 


Working fluids such as ateam, air, and exhaust gases have 
considerable vibrational heat capacity at high temperatures. If these 
gases have relaxation times comparable with or shorter than the in- 


tervals during which temperature changes occur in the gas, losses 





? 
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must be expected. Tor turbine worcing fluids, Kantrowitz!!) eati- 
mates that in sore cases the losses at high temperatures due to this 
adjustment las can be comparable with the losses due to skin fric- 
tion. Unfortunately, there exists 2 scarcity of data on the marni- 
tudes of the relaxation times of most gases and, without these values 
being known at verious temperature levels, their effects cannot be 
analyzed. 

Various experimental methods have been devised to measure 
relaxation tires of gases. Of these, only the two discussed here 
have been used with any degrec of success. A third method involv- 


ing a shock tube has been proposec, but has as yet not actually been 


used. 
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WT. EXPERIMENTAL METHODS 


A. THE SOUND VELOCITY METHOD!) 

The sound velocity method can be adapted for the determina- 
tion of the relaxation times of gases up to temperatures of the order 
of 2300°K. The two essentials are (a) the production of sound waves 
of a definite frequency and (b) the accurate measurement of temper- 
ature. 


3) » consists essentially 


The apparatus, originated by Pierce! 
of a carbon tube which can be heated to a high temperature and in 
which a stationary train of sound waves can be set up from a quartz 
crystal vibrating piezo-electrically. A piston traverses the tube. 
The position of the piston is read on an accurately graduated steel 
rule. This piston functions as a reflector te set up a standing wave 
system inthe tube. The amplitude of vioration of the crystal changes 
as the reflector is moved along the tube. Resonance is detected by 
observing the magnitude of the plate current of the oscillator which 
maintains the quartz crystal in vibration. The half wave-length is 
determined by recording the different positions of the reflector at 
which resonance occurs. The temperature of the gas in the tube is 
measured by sighting an optical pyrometer of the 'disappearing fila- 
ment" type on to the face of the movable reflector. An clectric fur- 
nace surrounds the carbon tube so that the temperature may be ad- 
justed to higher levels. Figure 1 shows the general layout of the ap- 
paratus. 

The velocity of sound, V, in the gas is computed from the 


measured half-wave length in the tube and from the known frequency 
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of vibration of tne crystal. The ratio of specific heats is given by 


the formula, 


The Vm /RT, 7 (3) 


where M is the molecular wdight, R is the gas constant, and T is 
the absolute temp erature. The correction necessary to the above 
formula for departure from the perfect gas laws'is negligible, since 
at the temperatures involved most gases are essentially perfect 


ases. From the observed data, an epparent ratio of specific heats, 


eS 


which we may call y', is computed from Eq. (3). Also from the 


relationship, 


Cy/R =1/Cs-1) , (4) 


an apparent specific heat at constant volume, Co. can be obtained. 
The crystal used for the production of sound waves, if cut 
properly, can be made to vibrate in either the flexural mode or the 
longitudinal mode by suitable electronic circuitry. The two fre- 
quencies produced may vary by as much as 20,000 cps or more. 
The apparent specific heat obtained at each of these frequencies 
may differ considerably. Kneser'4) has put forth the following ex- 
planation for this phenomenon. As a sound wave of a given frequency 
passes through a gas, a part of the vibrational energy fails to follow 
the acoustic cycle, and the gas departs from its equilibriu:: condi- 
tion. In fact, a part of the vibrational energy vanishes from the 


expression for the adiabatic elasticity, thus giving an increased 
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ek 
sound velocity. Kneser's theory led to the following forrnula for 


the apparent specific heat, derived from the velocity of sound data, 


iS %etwenreu 
vo —— oe | . as (5) 
Cy + 417 n T Ca 
where Cy, = true specific heat 


n = frequency of the sound wave 
Ca = specific heat of translation and rotation (= ae for 
a diatomic molecule). 
The values of the apparent specific heate raay be substituted in this 


equation to give two equations at each ternpérature for the determina- 


tion of T and Cy. 


B. IMPACT TUBE METHOD 

In 1942 Kantrowitz'!) proposed the impact tube method for 
finding relaxation times and rnade sufficient rneasurements to veri- 
fy that the experimental results checked with the theoretical calcu- 
lations that were available at the time. Grigrith! >) later extended 
this work and in 1950 reported on the results of his investigations, 
listing the relaxation times determined and expressing confidence 
in the method. 

In the impact tube method, gas at rest is allowed to expand 
adiabatically through a nozzle, and the flow thus produced is brought 
to rest again by compression at the nose of an impact tube (pitot 
tube). The total head defect in the flow, as a measure of the en- 
tropy increase, is directly associated with the heat-capacity lag of 
the gas and can be used to determine the relaxation time. The pro- 


ces6 is shown schematically in Figure 2. 


#6. 

The underlying principles involved may be summarized 
briefly as follows. The equations of momentum, continuity, and 
energy for a compressible fluid with possible separate transla- 
tional and vibrational temperatures, together with the standard 


definition for entropy, may be combined to establish the rela- 


tionships, 
| 1 1 
ds = ls. 5 2| dq (Cyn Tin) ’ (6) 
2 = 
s ( CLT ie + CpT +q /2) “{ (7) 


The notation Cy, refers to the specific heat of vibrational 
motion. We have then that C Sib JE... oT. > and if CE is the 
specific heat at constant pressure excluding contributions from vi- 
bration, we also have that Cy» Cy Can. 

(6) 


Using the theory of Landau and Teller’ ", the variation of 


CHib Trib with time may be expressed as 
Z (ce) = €.2eT os Vis ; (8) 


on the basis of considerations involving molecular collisions, the 


relaxation time is found to have the following functional dependence: 


~*/3 ~¥3 °*s Ys 2 
7 ee p. (vs) (m/T) exp | Cvs) (m/T) | . (9) 
The quantities y , 5, and m represent, respectively, vibrational 


frequency, molecular radius, and molecular mass. 


The departure from equilibrium between vibrational and 
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translational temperatures was expressed by Kantrowitz as 
6 = Gag(Tee- )- (10) 


In this experiment the changes in T are small compared to T itself, 
and q. (6) is redaced to 
dy xa edt t- Tas 2 e‘dt 


= 1 
— 4 es ae YC opt ° al 


~~. 
T is taken as the average of TT,-.. over the whole process. 


The total head defect Qh is measured as the pressure dif- 
ference p,- Pa > and since this difference is small coxnpared to 
Pp. » the differential change in entropy, ds = - Rdp/p, may be closely 


approximated by 


AS = AhR / po : (12) 


Combining Eqs. (11) and (12) we obtain the relationship, 


2 


dt 
Ah/p, = 2 : (13) 


If the temperature is eliminated from Eqs. (7), (8), and (10), 


it is found that 











2 
de a Cpe me Cee d(q/2) 
= - ; 14 
dt Cpt Cp dt vai 


For the fast stopping compression process at the nose of the impact 
tube, the second term of this expression becomes small in compari- 
son with the other terms, and the change in € becomes 
% 
Cs A q 
e 
ee 2 


Ae = ° (15) 





-§. 
Since at the beginning of compression @€ is zero and the velocity 


if gq, » Eq. (15) becomes 














2 
_ Coin % 
<_< -— (16) 
, @ 
With q = 0, Eq. (14) is reduced to 
de a Cpe 
— 3 - (17) 
dt Cp T 
For instantaneous compression the result is accordingly, 
oO 2 
2 : u 
Ahi. edt Ce Cvib S: 


7 —_—" Cyt aeRO T | 2Cp 


Aft this point it is convenient to introduce the following set 


of dimensionless parameters: 
(89s ‘ e x 3 
T = Cpeqt/Ced , € = 2CpE /Cyoa, , 


t mers ; q 4/4 ; 


where dis tne diameter of the impact tube. 


(19) 


lf Eqs. (13), (18), and (19) ave now combined, it is found 


that 
Ah 2 
a and 


Ah: — 


, ( , ‘ 2 « 
In order to determine the value of € in this expression it is neces- 


sary to solve Eq. (14) exactly, the integral being 


t, c, t, (3) 
e = exp] - Code exp Codt | Cyto DLA? ge 5 (22) 


which, stated in terms of the dimensionless parameters, icads to 


the expression 





—= 


dt’ dq dt. (32) 
dt! 


The integration of Eq. (22) has been done numerically for 
various values of T'. The integration of Lq. (20) has been car- 
ried out graphically, and from these resulting values a curve of 
ah/ Ah; vs. T' has been prepared. With this curve on hand an 
experiment may be conducted using any gas, and the relaxation 
time may be determined from the measured values of TY, be 9 
and p, . The vaines of ene F Cp ; Cp , and R for thé gas can 
be found by the use of well-imown statistical raecthods. 

The procedure then is as follows. 

a) Calculate ‘T’ from the relation ship, T= Te (Cp. / pe) 

b) Calculate q, from the relationship, Q, = [2 Cp CT -7,)| 

c) Cdiciiste Oh, from isq. (28). 

adj) Caleulate ah from the reiationship, Ah = Po -P3. 

e) Use the ratio of ah / Ah; to enter into the graph of 

ah/aoh,. vs. T to find the corresponding value for wpe 


f£) Calculate the relaxation time, T , from the relationship, 


v= ¥ Cpd/Cpq,. 


C. EXPERIMENTAL RESULT 


Table I liste the values of the relaxation times of twenty cases 


determined experimentally by various investigators. The informa- 


tion in this table was obtained, for the most part, from Reference (5). 
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WI. THEORETICAL CALCULATIONS 


The accumulation of a large body of experimental data led to 
efforts on the part of investigators to determine vibrational relaxa- 
tion times using theoretical calculations. The method of Landau and 
Teuier”), aeveloped in 1936, and xelined by Bewhe and Tenes! Gp, 
1940, has received general acceptance. A more detailed method de- 
veloped by Schwartz, Slawsky, and Herzfeial® 2) was published in 
1952. This method of calculation is based on the fact that an equiiib- 
rium between the various degrees of freedom can be reached only 
by means of a transfer of energy taking place during a molecular col- 


lision. Thus for a pure gas it can be said that 
 =M,,P 
= aa : (23) 


where M a8 is the number of collisions a molecule experiences per 
second and is proportional to the density, while P is the probability 
that upon a collision a quantum of vibrational energy will be converted 
to translational energy. For a gas composed of polyatomic molecules 
or mistures of dittcrent types of molecules, there arises the possibility 
of an energy exchange between the different vibrational degrees of 
frecdom during a collision as well as exchanges with the translational 
degrees of freedom, and a separate relaxation time Paiiet be ascribed 
to each of these energy transfers. 

The theory of Schwartz, Slawsky, and Herzfeld will be dene in 
this paper to compute the variation of vibrational excitation probabi- 
lities with temperature of mixtures of Os and No, and of HF and ti,. 


For a mixture of molecules of type a and b, following the notation of 





sie 
Schwartz et al, we may define the probabilities FP 7 (4,0) ; 
po (hb), FP, (a,b), Pha), (a,b), and (ba). 


> ovo se 


In each of these cases, the molecule, the symbol of which appears 
first in the parenthesis, was initially in the first excited vibrational 
state and finally in the ground ie zeroth vibrational state. The sub- 
script designates the behavior of the vibrational quantum number of 
the molecule, the symbol of which appears second in the parenthesis. 
Each of these probabilities has been tabulated and graphed as func- 
tions of ternperature for both the Os - Ny and HY - Hi, gaseous mix- 
tures. 

The probability a =. (b,3) may be obtained from the rela- 


tion 


-h(v,-P,) JKT 
a (b, 3) = P. (4,6) e of . (24) 


The others are determined independently. These probabilities may 


then be used to obtain the relaxation times from the relations 


= . 46) \1i- - 5 
Tae) ae ery ) (25) 
and 
—. mM (a,b) (1 - Site a ) 
=, e e 
Tas) a phe o-Fo (26) 


Here M, A and Mab denote the number of collisions of the designated 


type per second. We may define similarly VT {b,b) and VT (b,2). 


Say 


For the complex excitation process we may write 


MA 


4 et 
= XX, Mg, ,R.,G.6 (27) 








a= 
with a similar definition for Tg (h.4)} Here KMyand Ay are respec- 
tively the mole fractions of species a and species b in the mixture. 
It should be noted that only under special Chacumnetant ak can the re- 
laxation times defined above be grouped together to represent an ef- 


fective relaxation time which corresponds to a measured relaxation 





time. 
In general, we may write for the probabilitice 
= | 
P= CP , 128) 
where 
SuTPAE 2 Vo -S 2.. 
Dp! => 4.13 a 0, e {|v G) 5 (29) 
a” hh ‘ 
" “Wa -E/«T 
Cz 045 | +e") a a| e ; (30) 
“fe 2. 
o = Sp¥s : (31) 
2kT 
sr = Spe” + AE {32) 
2kT 2kvT ° 
=! 
v, = autkTac \S_ ae swaTar) (33) 
° a a 
a” hp 2p oh a 
and 
— . . ES 
A a (34) 


Here » is the reduced mass of the two colliding molecules, AE 
represents the amount of energy exchanged between vibration and 
transiation, ¥% and € are the usual constants in the Lennard-Jones 


interaction potential, and Vi) is the customary perturbation 
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integral over harmonic oscillator wave functions. For molecules 
whose value of € is not too great, it is aufficient to write 
o* = 17.5 /ve . In the case of HF, where sufficient virial 
coefficient, or transport data are not available, the collision di- 
ameter, WY , was estimated by plotting re as a function of the 
electronegativity of the halogen atoms. From Fig. 3, it can be 
seen that the value of Ve for HEF is approximately 1.7A. 


The condition 
M a *fe > Se (35) 


will not hold below a certain temperature for each of the probabili- 


ties. ‘We need this inequality to write the expansion 


2. 
i ior = pav3’ 2pve® iy : hs 


It is not necessary that 
pvs*/2 >> AE 


although the larger the inequality, the better the convergence. 
Schwartz and Herzfeld, in the appendix of their paper, have given 
amore complex method of computation which applies in the region 
of or below the limit ternperatures defined by Eq. (35). 

The following derivation applies for the determination of the 
temperature limit for which the basic equations used here wili be 


applicable. The condition required by inequality (35) may be re- 


‘eo > (20e/p)”™. 


If we combine this expression with Eq. (33) we obtain 


stated as 
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ant l3 =" fe 
(és KTA€ | _ AE. AE (ax KTAE Ey" > aaa 
- Ta (37) 








o* hp 2p \ otha 
H 2. 
4m kAE _ 
oh pA 3 eb -am Aaey 6 : 


inequality (37) is reduced to 
Ys nt /s 
(IT)? -e/4(IT)® > B , 


whence 
yar - er - B/43% >o., (38) 


We may obtain from this inequality the lower temperature limit of 


‘TT such that 
3 
Vises 1.758 6 3 ; 
124 wh /ae\*® 
Tiimit $ Tt (=) ° (39) 


Eq. (39) has been used to compute the values of the limit tempera- 
tures shown in Table I. 

A compilation of the constants used in making the calcula- 
tions indicated by the basic equations is shown in Table HI. The 
constants in the Lennard-Jones interaction potential, 7 and % 
were obtained from data listed in the ''Transactions of the ASMiz (19) 
with the exception of the value for % for HF which was es  atadehi 
in the manner already described. The vibrational constant V was 
found in each case from ''Spectra of Diatomic Molecules" by Herz- 


berg’ !?), 
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The results of the indicated computations are listed in Table 
IV and shown graphically in Figures 4 and 5. The valuee indicated 
for the O, ~ N, mixture agree within a few percent with those listed 
in the paper by Slawsky and Herzfeld. It is believed that where dif- 
ferences do exist, they may be attributed to a variation in the selec- 
tion of force constants and low-velocity collision diameters inagmuch 
as the published values are not always consistent. The calculated re- 
laxation tirnes also agree reasonably well with values determined 
experimentally, as discussed in the paper of Schwartz et al. No 
data are available for comparison with our results for the HF - He 


probabilities. 
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IV. VIBRATIONAL LAG IN THE HYDROGEN-F LUORINE 


ROCKET MOTOR 


It is now possible to present some bricf considerations in 
determining whether or not the adiabatic expansion through the 
de Laval nozzle in the Hy - F, rocket motor may be taken as vi- 
brational near-equilibrium flow. The reactants are taken in such 
proportion that Hs will be in excess, so that the products of com- 
bustion passing through the nozzle will be mostly a mixture of 
HF and H > 


Following a procedure analogous to that of Penner!! 2), in 
his treatrnent of chemical reaction in nozzle flow, we may take as 
an estimate for the vibrational temperature lag for near-equilibrium 


flow that 


AT = (oT/ot)r. ris 


For representative nozzles in 1,000 psia thrust chambers, it has 
been found!) that the cooling rate (-DT/Dt) = (Tc-Te)fey* 3x107“ 
°K/sec. The use of very small nozzles will, of course, increase 
the cooling rate above this valuc; however, for nozzles in larger 


-/(/ O-@ 
i</sec, 


thrust engines, (-DT/Dt) will generally be less than 3 x 10 
thereby leading to better approximations of near-equilibrium flow 
than we estimate below. 

From Figure 5 it +e be seen that the probability P _,, (a,b) 
is approximately as large as some of the other probabilities and 


much larger than others. The smaller probabilities contribute 


little to the effective relaxation time and may be neglected. The 
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probabilities, which are approximately as large as P ,,G,pb) ; 
all contribute to the effective relaxation time. For the purpose of 
Simplicity, however, it may be assurnec that the inaximum possible 


value of the effective relaxation time can be approximated as 
Tw, = Tap & 1/X,X,Mg,R.,Ge) . (41) 


The actual vaiue of Y must, of necessity, Se smaller than TY as 

the other probabilities, which are approximately as large as = »(4,b) 

have not been considered in computing Y%, and wouid lower its 

value. li we can show, however, that vibrational near-equilibrium 

flow exists for the calculated upper limit, Ty » the conclusion 

will certainly hold for the real Y whichis smaller than Ty . 
Since Mab ,» which denotes the number of collisions per 

second between HF and Hs molecules, increases markedly with 

pressure, it is apparent that if the conditions at the nozzle esit 

are such that vibrational near-equilibrium flow exists there, then 

it must also exist everywhere in the nozzle. The value of AT 

was therefore computed at the nozzle exit, where it is a maximum. i 
In the present computations, the nozzic exit pressure was 

taken as one atmosphere. The product Xs x b will generally lie 

between the comparatively narrow limits of 0.25 and 0.20, and 

will usually be about 0.23. The product of Myy and P (ab) 

will not vary appreciably from 4x 10° between the a tievades 

1400 to 3000° K. Therefore, for sake temperatures within this 


“1 
range, Tr $ | 2x10) sec and, for (-DT/Dt) $3 x 10! °x /sec, 


*For amore cetailed description of iterative lag calculations dur- 
ing nozzle flow, see Reference (14). 
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AT, < 3°K, 
The estimate that over the length of the nozzle the vibra- 
tional temperature lag cannot exceed 3° K indicates that vibra- 
tional near-equilibrium flow occurs during expansion in a de 
Laval nozzle for representative hydrogen-fluorine rocket mo- 


tors. 
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A. Sonic V elocity Method 


Hydrogen 

Oxy gen 

Carbon monoxide 
Chlorine 

Carbon dioxide 
Sulfur dioxide 
Amrnonia 
nthylene 
Methane 

Carbon disulfide 
Nitrous oxide 


B. Impact Tube Method 


Hydrogen 
Nitrogen 
Nitrous oxide 
Carbon dioxide (comin) 
Sulfur dioxide 
Ammonia 
Mythl chloride 
Sthane 
Ethylene 
-ropylene 
Propane 
Ethylene oxide 
Eutane 
Butadiene 
Methane 
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287 
289 
291 
291 
293 
290 
295 
287 
289 
294 
294 
303 
294 
294 
289 


| 
| 
| 
| 
| 
! 





0.021 
0.602 
l.l2 
Ze OZ 
0.37 
G.i2 
G. 202 
0.003 
0. 207 
0.096 
0. 007 
1.23 
0.018 
3.013 
0.42 


ere: emer eee a 


TABLE Il. ixperimental results reported in Reierence 5 


for the relaxation times of gases 
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TABLE Il. Lower temperature limits for which 
the probabilities can be calculated from the theory 


of Schwartz, Slawsky, and Herzfeld. 
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TABLE I. Values of constants used in computations for 


the vibrational excitation probabilities. 
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rig, 1. Schematic drawing of apparatus for determining vibrational 


relaxation time by the sound velocity metnod. 
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Fig. 2. Schematic drawing of the experiment and temper=ture 
history of tne gas for determining vibrational relaxation time 


by the impact tube method. 
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FIGURE 3- COLLISION DIAMETER AS A FUNCTION 
OF ELECTRONEGATIVITY. 
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FIGURE 4— VARIATION OF (a=0,,b=N,) PROBABILITIES WITH 
TEMPERATURE. LIMIT TEMPERATURES CORRESPOND TO 
THE JUNCTION BETWEEN THE SOLID PORTION AND THE 
DOTTED PORTION OF THE CURVES. 
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FIGURE S- VARIATION OF (o=HF,b=H,) PROBABILITIES WITH 
TEMPERATURE. LIMIT TEMPERATURES CORRESPOND TO 
THE JUNCTION BETWEEN THE SOLID PORTION AND THE 
DOTTED PORTION OF THE CURVES. 
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